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ABSTRACT 

The large relative sizes of circumstellar and circumplanetary disks imply that they might be seen in 
eclipse in stellar light curves. We estimate that a survey of '^10'* young (^^10 million year old) post- 
accretion pre-main sequence stars monitored for ~10 years should yield at least a few deep eclipses from 
circumplanetary disks and disks surrounding low mass companion stars. We present photometric and 
spectroscopic data for a pre-main sequence K5 star (ISWASP J140747.93-394542.6 = ASAS J140748- 
3945.7), a newly discovered ~0.9 M© member of the ~16 Myr-old Upper Centaurus-Lupus subgroup 
of Sco-Cen at a kinematic distance of 128 ± 13 pc. This star exhibited a remarkably long, deep, and 
complex eclipse event centered on 29 April 2007 (as discovered in SuperWASP photometry, and with 
portions of the dimming confirmed by ASAS data). At least 5 multi-day dimming events of >0.5 mag 
are identified, with a >3.3 mag deep eclipse bracketed by two pairs of '^1 mag eclipses symmetrically 
occurring ±12 days and ±26 days before and after. Hence, significant dimming of the star was taking 
place on and off over at least a '^54 day period in 2007, and a strong >1 mag dimming event occurring 
over a ~12 day span. We place a firm lower limit on the period of 850 days (i.e. the orbital radius 
of the eclipser must be >1.7 AU and orbital velocity must be <22 km/s). The shape of the light 
curve is similar to the lop-sided eclipses of the Be star EE Cep. We suspect that this new star is 
being eclipsed by a low-mass object orbited by a dense inner disk, further girded by at least 3 dusty 
rings of optical depths near unity. Between these rings are at least two annuli of near-zero optical 
depth (i.e. gaps), possibly cleared out by planets or moons, depending on the nature of the secondary. 
For possible periods in the range 2.33-200 yr, the estimated total ring mass is ~8-0.4 M-Mooni^i the 
rings have optical opacity similar to Saturn's rings), and the edge of the outermost detected ring 
has orbital radius ^0.4-0.09 AU. In the new era of time-domain astronomy opened by surveys like 
SuperWASP, ASAS, etc., and soon to be revolutionized by LSST, discovering and characterizing 
eclipses by circumplanetary and circumsecondary disks will provide us with observational constraints 
on the conditions which spawn satellite systems around gas giant planets and planetary systems 
around stars. 

Subject headings: binaries: eclipsing — planets and satellites; formation — planets and satellites: 
rings — stars: individual (ISWASP J140747.93-394542.6, ASAS J140748-3945.7) 
stars: planetary systems — stars: pre-main sequence 



1. INTRODUCTION 

The radii of circumstellar and circumplanetary disks 
can vastly exceed those of stars and planets. A com- 
panion star of a young stellar binary system can host 
a circumstellar disk, and likewise a giant planet in a 
young stellar system can host a circumplanetary disk. 
Because the disks are large, the probability that a ran- 
domly oriented system exhibits eclipses may not be neg- 
ligible. These disks are particularly interesting as they 
could be seen in eclipse during the epoch of planet for- 
mation (in the ca se of a companion circumstellar disk; 
IGalan et al.l l2010l ) or during the epoch of satellite for- 
mation (in the case of a circumplanetary disk). In this 
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paper we consider the possibility of discovering eclipses 
by dust disks of low-mass companions in long period or- 
bits. With the advent of long term and large scale pho- 
tometric surveys, strategies can be developed to discover 
young systems eclipsed by disks. 

Some well-known long period eclipsing systems have 
been interpreted in terms of occulting dark disks as- 
sociated with an orbiting c ompanion, w ith the best 
examples being e Aurig a e (iGuinan & DcWarf 200^; 
iKloDDenborE et al.' '2010; Chadima et al.l l2011l\ EE 
Cep ( Mikolaiewski &: Graczyk 1999: G r a.czvk et al.ll2003l : 
'Mik olaiewski et all r2605i: iGalan et all I2010D. aiid the 
newly identified OGLE-LMC-ECL- 17782 (jGraczvk etahl 
1201 ll i. EE Cep exhibits long (30-90 day) asymmetric 
echpses with a period of 5.6 years and depth of ~0.6 - 2.1 
magnitudes. The primary obje ct is a B5e giant star and 
only primary eclipses are seen ("Mikolaie wski fc GraczvlJ 
,1999 ). Structure seen in the wings of the eclipse has 
recently been interpreted in terms of rings and gaps in 
a formi ng planetary syste m around a lower mass sec- 
ondary (jCalan et al.ll20ldl ). e Aurigae is the eclipsing 



2 



Mamajek et al. 



system with the longest known period of 27.1 years. Only 
primary eclipses are seen and they last almost 2 years. 
While the mass of the object hosting the dark occult- 
ing disk exceeds that of visible F star, the masses of 
the two stars are not well constrained. However, emis- 
sion lines and UV emission suggest that the hidden ob- 
ject i s a B-type star (see discussion by iChadima et al.l 
I20T1I) . Infrared emission from the disk was detected 
with IRAS (IBackman fc GillettI ITOSSl) . OGLE-LMC- 
ECL- 17782 (MACHO J053036.7-690625; 'Graczv k et all 
|2011| ). is a 13-day eclipsing binary in the Large Mag- 
ellanic Cloud that demonstrates wide, flat-bottomed 
eclipses like e Aur, but there are changes in the light 
curve from eclipse to eclipse, and tr a nsieii t features vis- 
ible at other phases. iGraczvk et al.l (|2011[ ) suggest that 
the system is a detached binarj0 where the secondary is 
"partially hidden within a semi-transparent, dark, elon- 
gated body or disK^ and there are likely "transient struc- 
tures in the system (disk debris?) responsible for addi- 
tional minima at different orbital phases when one of the 
stars is hidden behind them.'''' 

In this paper we present the discovery of a solar mass 
pre-main sequence non-accreting star exhibiting a long, 
unusual eclipse similar to those seen in EE Cep and e Au- 
rigae. The mass of the star and lack of detected infrared 
emission suggest that the host object for the eclipsing 
disk is low mass, possibly a low mass star, brown dwarf, 
or giant planet. Hence we consider both circumsecondary 
and circumplanetary disks as possible occulting objects. 

Why should one also consider circumplanetary disks 
associated with giant planets? First, natural satellites 
are a ubiquitous feature among the giant planets in our 
solar system, and most likely among extrasolar gas gi- 
ants. The existence of such satellites, and the H/He-rich 
atmospheres of gas giants hint these planets likely formed 
with large gas and dust disks {rpianet « rdisk < rnui) 
that were originally accretion disks feeding from circum- 
stellar material. These disks would then evolve passively 
after the circumstellar reservoir was depleted (Sec. 2), 
with matter accreting on to the planet, grain growth and 
proto-satellite accretion, and depletion through other 
mechanisms (e.g. Foynting-Robertson drag, radiation 
pressure, photoionization, etc.). 

A simple thought experiment illustrates the poten- 
tial observability of moon-forming circumplanetary disks 
around young gas giants (and indeed this was the back of 
the envelope calculation that spawned our interest in the 
interpretation of the eclipsing star discussed in Sec. 3) . If 
one were take the Galilean satellites of Jupiter and grind 
them up into dust grains, and spread the grains uni- 
formly between Jupiter and Callisto's orbit, one would 
have a dusty disk of optical depth O(IO^). The size of 
such a proto-moon disk in this case would be a few solar 
radii - i.e. large enough and optically thick enough to 
potential eclipse a star's light. Of course such a disk need 
not be face on - more likely the disk would have a non- 
zero inclination respect to the planet-star orbital plane, 
so the star need not be completely geometrically eclipsed 

5 The author d ereddened the UBV photometry for OGLE-LMC- 
ECL-17782 from [Massed pOOl . catalogued as M2002 #148104) 
using the Q-method, and finds that the primary is most likely a 
slightly evolved ~B2 star with (B-V)o ~ -0.23 and My ~ -3.7. 
However this calculation does not take into account the difference 
in metallicities between LMC and local Galactic massive stars. 



by such a circumplanetary disk. The rings of Saturn have 
optical depth near ~1 even at a relatively old age (4.6 
Gyr), however the vast majority of mass orbiting Sat- 
urn is locked up in satellites {Mrings — 10"'* M satellites)- 
Fresumably a disk of much higher optical depth and sig- 
nificant radial substructure existed during the epoch of 
satellite formation. While there have been studies in- 
vestigating the detectability of thin, discrete planetary 
rings similar to Saturn's (e.g. iBarnes fc FortnevI [200I 
lOhta et"al] '20091 , there has been negligible investigation 
of the observability of the dense proto-satellite disks that 
likely existed during the first ~10^ years. Relaxing the 
assumptions about the size, mass, composition and struc- 
ture of the disk in our back-of-the-envelope calculation 
has little impact on the feasibility of the idea that dusty 
disks of high optical depth may be a common feature of 
young gas giant planets, and such objects may be observ- 
able via deep eclipses of young stars. 

We first estimate the timescale of an eclipse by a cir- 
cumsecondary or circumplanetary disk in Sec. 2. This 
is done first so that the photometric light curve of our 
object can be interpreted in terms of eclipse models. We 
then present and discuss the properties of our candidate 
long period eclipsing system in Sec. 3. In Sec. 4 we dis- 
cuss the probability of detecting eclipses using time series 
photometry of large samples of young stars. A discussion 
and summary follows in Sec. 5. 

2. CIRCUMSECONDARY AND CIRCUMPLANETARY DISK 
ECLIPSES 

The multiplicity of class I young stellar objects (YSOs) 
in embedded clusters and pre-MS stars in young associ- 
ations is high and ranges from 20 to 60% for the ~10*'^- 
10^'^ AU separation range and mass ratio s of '-^O.l-l (see 
the comprehensive rev iew by [Puchcnc e t al.l 120071 and 
also iKraus et al.l 120 lit ) . Among the binaries found in 
young clusters , stars in diffe rent stages o f disk evolution 
are n o t rare (iHartig an & KenvonI l2003t iMcCabe et al.l 
I200I iMonin et ahf l2007t IFrato fc Weinbergeii l2oTor 
These binaries, known as "mixed pairs", have one star 
hosting a disk or actively accreting and the other lacking 
a disk or signatures of accretion. While the two most well 
known long period eclipsing disk systems (EE Cep and 
e Aurigae) involve massive stars, the large fraction of bi- 
nary systems in young associations and clusters suggest 
that long period eclipsing circumsecondary disk systems 
may be discovered in lower mass systems. 

Recent explorations of circumplan etary disks separate 
the disk evolution int o two phases ([Alibert et al.ll200"5l : 
iWard fc GamiDll2010D . In the first phase the circum- 
planetary disk is fed by material from the circumstellar 
disk and the circumplanetary disk acts like an accretion 
disk. In the second phase the circumstellar disk has dis- 
sipated, and the viscosity of the circumplanetary disk 
drives both accretion onto the planet and causes the disk 
to sp read outwards. The lack of differentiation of Cal- 
hsto (jAnderson et alll200l|) suggests that the accretion or 
formati on times cale for all the Galilean satellites was pro- 
longed (jCanup fc Ward 2002) and would have happened 
during the second p hase of evolution after the circumstel- 
lar disk dissipated (iCanup fc Wardll2002l: lAlibert et al.l 
I2OO5I IWard fc Canudl2010l:IMosaueira' et al JI2O1O0. Also, 
modeling of lapetus, the outermost regular satellite of 
Saturn, suggests that it formed 3-5 Myr after the for- 
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mation of Ca-Al inclusions ()Castillo-Rogez et al.ll2009t ). 
Since lapetus survived Type I migration it must have 
formed near the end of substantial accr etion onto Sat- 
urn from the circumsolar nebula. Mosquc ira et alJ ()2010l ) 
suggest that lapetus' large separation from Saturn's 
principal satellite Titan is suggestive that the remnant 
circum-Saturnian nebula may have had two components: 
a dense inner disk that spawned most of Saturn's regular 
satellites, out to the centrifugal radius near Titan, and a 
disk of much lower density beyond the centrifugal radius 
out - perhaps to Phoebe's orbit. Because of the extended 
estimated circumplanetary disk lifetime, we can consider 
the possibility that circumplanetary disks can be seen in 
eclipse against a young central star after the dissipation 
of the circumstellar disk. 

We consider two bodies with masses mi , rn2 in a cir- 
cular orbit with semi-major axis as and mass ratio 
In the case of a stellar binary the more 



u = — — — 

^ m\-\-m2 ' 

massive star is mi and the secondary is m2. In the case 
of a system with a single planet, m2 is the mass of the 
planet and mi the mass of the central star. The masses 
and semi-major axis set the Hill or tidal radius of the 

secondary = as (f )^^^- 

A disk surrounding m2 is described with two param- 
eters, the radius at which its optical depth is of order 
unity, r^, and the obliquity or axial tilt of the disk sys- 
tem, e with respect to the axis defining the orbital plane. 
The angle e is zero for a disk that lies in the orbital 
plane. It is convenient to define a size ratio ^ = rd/rn 
that represents the size of the disk in Hill or tidal radii. 

Studies of giant planets have made estimates for the 
size ratio £. Based on a centripetal radius argument 
IQuillen fc Trilling! (|1998[) estimated that an accreting cir- 
cumplanetary disk would have ^ = 1/3. Hydrodynamic 
simulations of planets e mbedded in circumstellar disks 
also can find ^ ~ 0.3 (jAvliffe fc Ba te 2009 ). A tidal 
trunc ation argument suggests ^ ~ 0.4 ()Martin fc Lubowl 
120 lit) . Theoretical models for Jupiter's circumplane- 
tary disk accounting for Galilean satellite formation af- 
ter the dissipation of the circui nstellar disk estimate 
smaller radii of £ ^ . 1 - 0.2 (ICanun fc WardI [20?i2l: 
IMaeni fc Coradia[200l IWard fc Canudl2010[) . For ref- 
erence, the outermost of the Galilean satellites, Callisto, 
currently has a semi-major axis that is only a small frac- 
tion of Jupiter's Hill radius a k, 0.0355ri/. Jupiter's Hill 
radius is about 743i?j where Rj is the radius of Jupiter. 
Estimates of the Hill radii th and the ratio of the orbital 
radii to Hill radii {rgat/TH) for the outermost large, regu- 
lar satellites for the giant planets in our solar system are 
summarized in Table [T] Column 6 of Table [T] estimates 
the timescale for a circumplanetary disk with outer ra- 
dius equal to the orbital radius of the outermost regular 
satellite to transit in front of the Sun, given the planet's 
mean orbital velocity. 

The results suggest that substantial circumplanetary 
disks of sufficient surface density for satellites to accrete 
must have existed around the giant planets in our solar 
system for some period during the post-T Tauri phase for 
our Sun. These circumplanetary disks likely had outer 
radii of £ > 0.01-0.05, and could have been detectable 
by eclipses with ^1-10 day timescales to observers along 
opportune lines-of-sight. 

A disk in a binary system might extend all the 



TABLE 1 

Hill Radii and Orbital Radii of Outermost Regular 
Satellites for Large Planets 
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(4) 


(5) 


(6) 


Planet 


Rh 


Vorb 


Outermost 


Tsat/rH 


^transit 




(AU) 


kms~^ 


Reg. Sat. 




(days) 


Jupiter 


0.36 


13.1 


Callisto 


0.035 


3.3 


Saturn 


0.44 


9.7 


lapetus 


0.054 


8.5 


Uranus 


0.47 


6.8 


Oberon 


0.008 


2.0 


Neptune 


0.78 


5.4 









Note. — Column (2) is the Hill radius th in AU, (3) is the 
planet's mean orbital velocity in kms~^, (4) is the name of the 
outermost large regular satellite for each planet, (5) Tsat/rH is 
the ratio of that moons orbital radius to the Hill radius, and (6) 
is the timescale ttransit that a circumplanetary disk with outer 
radius equal to the orbital radius of the outermost large regular 
satellite to transit in front of the Sun, given the planet's orbital 
velocity. Neptune does not have a system of large regular 
satellites. Its principal moon Triton is in a retrograde orbit, 
and w as likely captured as a c omponent of a binary dwarf 
planet IjAgnor fc Hamiltonll200g ^ 



way to its Roche radius (iLin fc Papaloizoul IT993I : 
lArtvmowicz fc Lubowl 119941 : [Andrews et al.l |2010D with 
^ ^ 1. Fo r example, the disk of HD141569A 
(jClampin et all [200l may extend to its Roche ra- 
dius when the second ary is at pericenter and ap- 
proaches HD14156 9A (|Augereau fc Papaloizoul [200l 
IQuillen et aL|[2005f ). The truncated disks of HD 98800 
and Hen 3-600 are consistent with tidal truncation 
(|Andrews et al.llMot) . 

A circular orbit for m2 would have a circumference of 
^ 2'KaB ■ The disk extends a distance ir^ sin e above 
and below the orbital plane. Thus the area of a cylinder 
that could be intersected by an eclipsing line of sight is 
A ^ An a Br d sin e. Here we have neglected the thickness 
of the disk at low obliquity. If the disk is seen edge on 
then A ^ inaphd where hd is the scale height of the disk 
near its opacity edge. We can define a factor 



y(e) — max I — ,sme 



(1) 



To estimate the probability that a system containing a 
disk is oriented so that it would exhibit eclipses we divide 
this area by Anap 



Porient 



(2) 



This probability is independent of the binary or 
planet's semi-major axis. 

Studies of primordial circumpla netary disks have cal- 
cula ted their thermal structure (|Canup fc Ward! 120021 : 
Alibe rt et al.l I2005D . The circumplanetary disk aspect 
ratio h/r is predicted to be in the range 0.1 - 0.3 where 
h is the vertical scale height a radia l distance r fro m the 
planet's center (e.g., see Figure 9 : 1 Alibert et all 12005). 
Hydrostatic equilibrium relates the temperature T of a 
gaseous disk to its vertical scale height and the speed 
D, of objects in orbit {h ^ Cs/fl, where sound speed Cs 
= (fcr)-'^/^(/xm//)~^/^), k is Boltzmann's constant, n is 
the mean molecular weight, and tuh is the mass of hy- 
drogen). Due to the low circular velocities for objects 
in orbit about a planet, a gaseous circumplanetary disk 
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would not have a small aspect ratio until its gas has dis- 
sipated. 

The temperature and scale height of a circumstellar 
disk is set from the dominant source of heat, which is 
from absorption of stellar radiation unless the accretion 
rate is high. A good rough guideline covering a wide 
range of conditi ons is aspect ratio /i/r ~ 0.05 - 0.1 (e.g., 
see Figure 1 by lEdgar et a"I]|2007[ ). The factor j/(e) also 
depends on the orientation of the disk. Strongly mis- 
aligned disks should be brought into rough alignment 
by shocks associated with tidal torques in only about 
20 binary orbital periods (Bat e et al.ll20Ciol ). For mod- 
est misalignment angles, misalignment can persist over 
considerably longer periods. Disks have been imaged in 
wide binary T Tauri stars that are misaligned wi th the 
binary orbital plane (e.g., HK Tau; Stapelfcldt et al.l 
Il998f ). If the system is not a binary but a triplet or 
a quadruple system then misalignment may be mor e 
common (as discuss e d by iPrato fc Weinbergeil l2010t ). 
iPrato fc Weinbergeil (|2010f ) emphasize (see their section 
4) that "even binaries with separations of a few tens of 
AU - or less - cannot be assumed to harbor aligned disks 
coplanar with binary orbits." Hence the obliquity distri- 
bution for circumsecondary disks may be wide. 

Evaluating eqn. [5] with illustrative parameters we esti- 
mate that the fraction of randomly oriented systems (at 
a given mass ratio) hosting a disk that are oriented so 
they could exhibit an eclipse is 



for^ent ~ 0.004 ( ^ 



1712 

mj 



mi + 1712 

Mr. 



y 

0.3 



.(3) 

Here we have taken y to be the mean of the distribution 
of y{e) that depends on the obliquity and disk height dis- 
tribution and Mj is the mass of Jupiter. Because of the 
large size of a circumplanetary or circumsecondary disk 
the fraction f orient of objects capable of giving eclipses 
is not low, and that such eclipses will be seen. 

The timescale for the eclipse to occur will depend on 
the angular rotation rate of the orbit 



eclipse ' 



G{mi + TO2) as 

\ — 



(4) 



' 1.5 days 



0.2 



mj 



VlAU/ 



where we have assumed that the radius of the occulted 
object i?i <g; rd and P is the orbital period of the planet 
or binary. 

The fraction of the orbit spent in eclipse is 



feclips 



eclipse 

p 



1/3^ 



-1/3 



(5) 



and is independent of the semi-major axis. 

Once a disk eclipse candidate is identified one could 
search for reflected star light from the disk. The area 
intersecting light from the star A 47rr^ siny(e)c where 
the order unity factor c depends on the orientation of the 
disk. The fraction of reflected star light would be 



(6) 



where the order unity factor c' depends on c, the disk's 
albedo and the dependence on scattering angle. The dif- 
ference in magnitude between the reflected light of the 
disk and 



5m - 9.8 - 2.5 log 



10 



0.2/ VlO-^ 



0.3 



•(7) 



This level of magnitude difference is not extremely high 
suggesting that it may be feasible to detect reflected light 
from circumplanetary disks with an adaptive optics sys- 
tem. 

2.1. Disk Lifetimes 

For young binary systems in nearby star-forming re- 
gions (with typical ages of <3 Myr) the estimated frac- 
tion of mixed systems with primary a weak lined T Tauri 
star and the secondary a classical T Tauri star is not 
low, and could be as large as 1/3 (jMonin et al.i r2007'). 
The mixed systems imply that the circumstellar disks 
around the primary and secondary can have lifetime s 
that differ by a factor of about 2 ()Monin et al.ll2007D . 
A tidally truncated disk around a low mass secondary 
is expected to have a shorter accretio n lifetime than the 
primary's disk (jArmitage et al.|[l99 9j). however there are 
other environmental factors such as dispersal of the host 
molecular cloud, birth cluster density, and disk evapo- 
ration that can influence multiplicity and disperse disks 
()Monin et al.l [20071 : IPrato fc Weinberg er 2010). 

We can estimate the lifetime of a circumplanetary disk 
by scaling from models for the proto- Jovian disk. The 
orbital period of a particle in orbit about a planet with 
radius near the Hill radius is approximately that of the 
planet in orbit about the star. The lifetime of a circum- 
planetary disk likely depends on the orbital period at its 
outer edge and so depends on f'^/^ times the planet's or- 
bital period. If ^ is similar for different circumplanetary 
disks, the lifetime of the secondary phase of these disk 
(after circumstellar disk dissipation) should be propor- 
tional to the planet's semi-major axis to the 3/2 power. 
We can use this scaling relation to estimate the lifetime 
of circumplanetary disks at larger distances from the 
star than Jupiter. The lifetime of the second phase of 
Jupiter's circu mstellar disk is estimated to be of order a 
million vears (|Canup fc Wardll200l lAlibert et al.l 120051 : 
IWard fc Canudl2010D . The lifetime of a circumplanetary 
disk at 25 AU we estimate could be about 10 times longer 
(or of order 10*" years) and at 100 AU (such as Fomalhaut 
b) a 100 times longer or of order 10^ years. 

The planet Fomalh aut b has been detected in two 
visible bands only (K alas et all I2008D . While the de- 
tected object is in orbit about Fomalhaut it's color is 
consistent with that of reflected light from the star. 
These observations led Kalas et al. to suggest that 
the p lanet hosts a circumplanet ary ring akin to Saturn's 
(see I Arnold fc Schneideiil2004[ ). which would extend to 
at least 20 Jupiter radii for an assumed albedo of 0.4 
to recover the observed fluxes. Another possibility is 
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that this is reflected light from a gaseous circumplan- 
etary disk. Due to the large semi-major axis of the 
planet, (and so large Hill radius) the lifetime of this disk 
would exceed that estimated for Jupiter's circumplane- 
tary disk, making this possibility more tractable. The 
light from Fomalhaut b is unresolved so the emitting ob- 
ject must be confined to a region smaller than the Hubble 
Space Telescope Advanced Camera point spread function 
PSF) full width at half- maximum (FWHM) of 0.5 AU 
Kennedv fc Wvattll20H"l ). The ratio of 0.5 AU to the 
planet's semi-major axis of 119 AU gives a constraint 
= 2.2 X 10"''. For ^ — 0.2 this gives a planet mass 
ratio of 2 X 10~^ and for ^ = 0.1 of 2 x 10"** so planet 
masses ranging fro m Saturn to N e ptune mass and lower 
than estimated by | Chiang et al. l ( |2009D but consistent 
with that predicted by QuiLleri (2006). 

3. A CANDIDATE ECLIPSING DISK 

3.1. The Star 

We are currently conducting a large scale spectroscopic 
survey for new low-mass members of the Sco-Cen OB 
association (Pecaut & Mamajek 2012, in prep.) using 
the RC spectrograph on the SMART^ 1.5-m telescope 
at Cerro Tololo. Sco-Cen is the nearest OB associa- 
tion to the Sun (mean subgr oup distances of d ~ 118- 
145 pc; Ide Zeeuw et al.l[T999f ) and consists of three sub- 
group s with ages of '^ll-l? Myr (T lPreibisch &: Mamaiekl 
|2008() . The survey sample consisted of ~350 stars with 
optical/near- IR colors consistent wit h having K/M spec- 
tral types, PPMX proper motions (jRoeser et al.l l2008[ l 
consistent with membership to the three Sco-Cen sub- 
groups, and X-ray emission det e cted i n the ROSAT All- 
Sky Survey (|Voges et al.l [1991 [2000l ). T he photomet- 
ric an d astrometric survey PPMX catalog (|R oeser et all 
120081) is complete down to V ~ 12.8 magnitude with 
typical astrometric accuracy of 2 mas yr~^. The sur- 
vey sample was cross-referenced with the stars with light 
curves in the the first public data release (DRl) of the Su- 
perWASP pubhc archivcQ (Butters et al.ll2010[ ). Super- 
WASP is a photometric sky survey for detecting transit- 
ing extrasolar planets with instruments in La Palma and 
in South Africa, which have continuously monitored the 
sky since 2004 a nd the DR l contains nearly 18 million 
light curves (PoU acco et all Hooei. Of our ~ 350 Sco-Cen 
survey stars, at least 200 appear to be new bona fide pre- 
MS stars, and SuperWASP light curves were available for 
138 of them. 

Among the SuperWASP DRl data for the new Sco- 
Cen members, wc|3 identified a star with a remarkable 
light curve (PPMX J140747.9-394542 GSC 7807-0004 
= ISWASP J140747.93-394542.6 = 3UC 101-141675 = 
2MASS J14074792-3945427; hereafter "J1407"). The 
light curve is dominated by a quasi-sinusoidal component 
with amplitude ~0.1 mag in the WASP-V photometric 
band, with periodicity of 3.21 days (consistent with rota- 
tional modulation of starspots, typical for young active 
stars), and a deep eclipse with maximum depth ^3 mag 
between HJD 2454213 (23 Apr 2007) and HJD 2454227 

® http:/ /www. astro. yale.edu/smarts/ 

http:/ /www. wasp. le.ac.uk/public/ 
* The deep SuperWASP eclipse for this star was first noticed 3 
December 2010 by Pecaut & Mamajek. 



TABLE 2 

Properties of Star 



(1) 

Property 


(2) 

Value 


(3) 
Ref 


«(J2000) 


14:07:47.93 


1 


<5(J2000) 


-39:45:42.7 


1 




-25.4 ± 1.4 masyr" 


-1 1 




-20.1 ±3.5 masyr" 


-1 1 


Spec, type 


K5 IV(e) Li 


2 


E(B-V) 


0.09 mag 


2 


Av 


0.32 mag 


2 


Dist 


128 ±13 pc 


2 


EW{Ha) 


0.2 A{emis.) 


2 


liW(Ll i AD7U7) 


0.4 A(abs.) 


2 




4500l^°° K 


2 


log{L/LQ) 


-0.47 ±0.11 dcx 


2 


R 


0.96 ±0.15 Rq 


2 


X-ray flux 


3.59 x 10-2 ct s" 


1 3 


HRl 


-0.04 ±0.42 


3 


log{Lyi/ Ltol) 


-3.4 


2, 3 




lO^^-S ergs-l 


2, 3 




3.20 days 


2 


Age 


~16 Myr 


2 


Mass 


0.9 Mq 


2 


Note. — References: (1) 'Zacharias et aL| 
pOToT). (21 this paper, C3)iVoKcs et al. (2000"). 

TABLE 3 
Photometry of Star 



(1) 


(2) 


(3) 


(4) 


Band 


Ao 


Mag 


Ref 


V 


0.55 fj.m 


12.31 ±0.03 


1 


I 


0.79 Aim 


10.92 ±0.03 


2 


J 


1.24 t^m 


9.997 ±0.022 


3 


H 


1.66 ^m 


9.425 ±0.023 


3 


K, 


2.16 ^m 


9.257 ±0.020 


3 


Wl 


3.4 ^m 


9.252 ±0.025 


4 


W2 


4.6 ^m 


9.276 ±0.021 


4 


W3 


12 fim 


9.141 ±0.033 


4 


W4 


22 /im 


8.907 ±0.388 


4 


Note. 


— (1) V-band is median of 



SuperWASP and ASAS measurements 
out of eclipse, given equal weight to 
both datasets, and the ±0.03 mag is 
a systematic uncertainty. The Su- 
perWASP photometry was conv erted 
to Johnson V using factors by Besselll 
pOOa) and assuming VsuperWASP = 
^Tych o {Fo llacco et al. 2006), ( 2) DE - 
NIS friie DENIS Consortium feosl') . 
(3) 2MASS (Cutri ct al. 20 01), (4) 
WISE first data release (, Wright et al.l 
[20Tch . 

(7 May 2007), with a complex pattern of roughly sym- 
metric dimming and brightening within ±26 days of 29 
Apr 2007 HJD 2454220. The properties of this star are 
listed in Table [5] and the relevant optical/IR photometry 
is listed in Table [31 We discuss the eclipse further in tj3.2l 
The spectral energy distribution for the star (plotted 
in Fig. [ij is consistent with a lightly reddened K5 star 
(E(B-V) = 0.09, kv ~ 0.32 mag), with no evidence 
for infrared excess ( via 2MASS and WISE preliminary 
release photometry; iSkrutskie et al.l 120061 : iWright et al.l 
l2010f ). Using a low-resolution red spectrum taken 
in July 2009 with the SMARTS 1.5-m telescope RC 
spectrograph, shown in Figure [2l we classify the star as 
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K5 IV (e) Li, i.e. a Li-rich K-star with neghgible Ha 
emission (0.2A equivalent width; i.e. "filled-in" ) , and 
Na doublet feature that is weaker than that for dwarfs, 
but not consistent with a giant eit her (so we adopt 
interm ediate luminosity class IV; iKeenan fc McNeill 
[l98l . 




Fig. 1. — Observed photometry for J 1407 taken from Table |3] 
(filled circles) compared to the spectral energy distribution for a 
lightly reddened K5 dwarf with E(B-V) = 0.09 (red dashed line). 
We have assumed that the Ks minus WISE band (Wl, W2, W3, 
W4) colors are zero. 



The star is located in the vicinity of the Upper 
Centaurus -Lupus (UCL) subgro up of the Sco 0B2 as- 
sociation (Ide Zeeuw et al.l[l999t ). and its proper motion 
is statistically consistent with moving towards the UCL 
convergent point (negligible peculiar velocity of 0.9 ± 1.8 
km/s), with a kinematic distance of 128 ±13 pc (simi- 
lar to other UCL membersj^. Using this distance, we 
place the star on the HR diagram, see Figure |3] (logT, 
logL/Lo = 3.66, -0.47). Factoring in the ±0.11 dex un- 
certainty in log(L/LQ), dominated by the distance un- 
certainty, the isochronal ages and their uncertainties are 
listed in Table ID Factoring in previou s age estimates 
for th e UCL subgroup (see summary in iMamajek et al.l 
l2002f). and new age estimates using th e F-star MS turn- 
on (jPecaut. Mamaiek. fc Bubar Il2011f ). we estimate the 
mean age of UCL to be 16 Myr with a ±2 Myr (68%CL) 
systematic uncertainty. The isochronal age for J1407 
is consistent with this value, hence we adopt the mean 
UCL age as the age for J1407. Three sets of evolutionary 
tracks predict similar masses for J1407: ~0.9 M© (listed 
in Table SD. 

The kinematics of the star, rotation period (as dis- 
cussed below). X-ray emission, and its preliminary HR 
diagram position, are mutually self-consistent with the 
interpretation that this star is a nearby (distance ^ 130 

^ The distanc e and peculiar velocity were calculated following 
[Mamaiek' (2005), using the UCAC3 proper motion for J1407, and 
the updated estimate of the mean space motion for UCL from 
[Chen et al.i (l2011l ): (U, V, W) = (-5.1±0.6, -19.7±0.4, -4.6±0.3) 
kms~^. 



TABLE 4 
Isochronal Age and Mass 
Estimates for J1407 



Age 


(2) 
Mass 


(3) 

Models 


Myr 


Mq 






0.90 ±0.08 


1 


97+15 


0.89 ±0.07 


2 


14lJ^ 


0.86 ±0.06 


3 



Note. — References for 
the models are a s f ollows : 

(1) IBaraffe et al.l JTglD, 

(2) Ec ss et al. (2000), Ji) 
ID'Antona fc Mazzitelli. (^1993). 



pc), ~ 10'' year old, solar- mass pre-main sequence (Pre- 
MS) star. 




5800 6000 6200 6400 6600 

A[A] 

Fig. 2. — Comparison of CTIO 1.5-m red optical spectrum 
of J1407 to CTIO spectra of four spectral standard stars from 
Keenan & McNeil (1989): TW PsA (K4V). N Vel (K5III), HD 
36003 (K5V), and GJ 529 (K6Va). 
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Fig. 3. — HR diagra m position fo r star J 1407 with isochrones in 
log(age/yr) from .Baraffe et al.l ( I199S ) overlain. 



3.2. Light Curves 
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Fig. 4. — Individual measurements of SuperWASP (top) and 
ASAS (bottom) V magnitudes for J1407 during early 2007. Ab- 
scissa is Heliocentric Julian Date minus 2450000. 1 Jan 2007 mid- 
night corresponds to HJD 2454101.5. The eclipse was seen in both 
photometric data sets. The eclipse was deep for about 14 days but 
is bookended by a gradual dimming covering a period of about ±54 
days. Long term median magnitudes outside of eclipse are plotted 
with dotted lines (V = 12.29 for ASAS, V = 12.45 for SuperWASP). 
The systematic difference is mostl y due to Super WASP- V being 
calibrated to the Tycho band iPollacc o et a l.| 12003), whereas 
the ASA S is converted to the Johnson V system via Hipparcos 
UPoiman ski 2002). 

The V-band light curve for J1407 during the year 
of 2007 from the All Sky Automated Survey (ASAS; 
(|Poimanskil[200l 'l and SuperWASP (Super Wide An- 
gle Search for Planets) surveys is shown in Figure |31 
The SuperWASP survey is an ultra-wide field (over 300 
sq. degrees) photometric survey is designed designed to 
monitor stars between V ~ 7-15 mag to search for 
transiting extrasolar planets (jPollacco et al.| [2006). The 
public data archive of SuperWASP photometry is de- 
scribed in Butters et al. (2010). The SuperWASP DRl 
photometry for J1407 contains photometry for approx- 
imately 29,000 epochs during 206 dates between HJD 
2453860 (2006.34) and HJD 2455399 (2010.56), with me- 
dian photometric precision of 0.023 mag. The light curve 
for J1407 from the SuperWASP datg^is dominated by 
(1) a sinusoidal component with amplitude ^ 0.1 mag- 
nitude in the WASP-V photometric band, with peri- 
odicity of 3.211 days (consistent with rotational mod- 
ulation of starspots, typical for young active stars, e.g. 
IMamaiek fc Hillenbrand 2008), and (2) a 14 day deep 
eclipse of depth > 3 mag between date HJD 2454213 (23 
Apr 2007) and HJD 2454227 (7 May 2007), bookended 
by a gradual dimmings and brightenings to the median 
brightness (Fig. |4]). The same light curve is also shown 
in Figure [5] using median nightly SuperWASP values. In 
Fig. ini we plot the SuperWASP light curve for a com- 
parison field star situated ^100" away from J1407, and 
of similar brightness (plotted over the same time period 
as J1407's eclipse and the same magnitude scale as in 
Fig. [5]). There is no evidence for similar complex behav- 
ior during the epoch of J1407's eclipse in the light curve 



for the comparison star. We discuss various scenarios for 
explaining the dimming of J1407 in 
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Fig. 5. — SuperWASP and ASAS photometry for the April-May 
2007 eclipse event (s). ASAS photometry is a single measurement 
each night, whereas the SuperWASP magnitudes are nightly me- 
dian values (the standard error of the median is plotted). The 
four (two pairs) peripheral dips are labeled and matched to their 
partner. The midway date for both the "A" and "B" dips coincide 
within a day of HJD 2454220 (29 Apr 2007). Dips Ai and A2 are 
'^51. 5 days apart, and dips Bi and B2 are ~24 days apart. Periods 
of low extinction are labelled Zi and Z2 and Yi and Y2 . 

The ASAS-3 archive also contained an extremely long 
time-baseline light curve for J1407 (ASAS J140748- 
3945.7), with photometry provided over 599 dates be- 
tween Feb 2001 and Sep 2009 (more specifically, HJD 
2451887 and HJD 2455088). The ASAS light curve is 
plotted with the SuperWASP photometry during the 
eclipse in Fig. 01 and the entire ASAS light curve for 
2001-2009 is shown in Fig. [71 ASAS data also shows a 
minor, but sustained dip in magnitudes between 2001.20 
and 2001.24 (~14 days) of approximate depth ~ 0.2 mag- 
nitudes (see Figure jT]). This magnitude difference is only 
2(7 above the night to night dispersion however because 
the dimming lasted a couple of weeks, the event stood 
out as unusual. If it was a true secondary eclipse, then 
the period should be 12.24 years and the next secondary 
eclipse would take place around 2013.46. High cadence 
photometry of J1407 in mid-2013 should be able to test 
the idea that the 2001 event might have been a secondary 
eclipse. 

A series of nightly V-band images were taken of J1407 
with the CTIO 1.3-m telescope in queue mode during the 
first half of 2011. Three consecutive 10 sec images were 
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Fig. 6. — SuperWASP light curve for the comparison star 
ISWASP J1407252. 03-394415.1 (GSC 07807-00572), a star of sim- 
ilar brightness (V = 12.24) to J1407 and situated 99". 6 away from 
it. The time span covers the same range as the light curve for J1407 
in Fig. \5\ 7829 photometric data points with median photomet- 
ric error ±0.018 mag are shown and 68% (95%) are within ±0.021 
(0.077) mag of V = 12.245. Despite the appearance of some discre- 

„ r\ent photometric points (with correspondingly large photometric 

Can be retrieved from http : //m . wasp . le . ac. uk/publlc/lc/in de^flVy„^ there is no evidence for any complex behavior similar to 
with the identifier ISWASP J 140747. 93-394542. 6 that seen for J1407. 
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Fig. 7.— ASAS light curve for J1407 between Feb 2001 and 
Sept 2009. An unconfirmed shallower eclipse might have occurred 
in early 2001. Magnitudes are reported for a given night, and 
connected by a thin dashed line 

taken nightly during 106 nights between 7 Feb 2011 and 
22 Jun 2011. Visual examination of the data, and com- 
parison of the brightness of J1407 to neighboring stars 
of similar brightness shows no evidence of for deep (>0.5 
mag) eclipses during this period. 

3.3. Eclipse Substructure 

While the deepest part of the eclipse is not well sam- 
pled in Figures|lland[5lthe eclip se of J1407 is asymmetric. 
iMikolaiewski fc Graczvkl (|1999( ) proposed that the asym- 
metry of EE Cep's eclipses were due to the disk impact 
parameter with the line of sight, and given the similari- 
ties of the central parts of their eclipses, we suspect the 
same for J1407. Nightly averages of SuperWASP data 
(see Figure [5]) exhibit two pairs of multi-day dips, la- 
belled Ai and A2, separated by ~24 days and Bi and 
B2 separated by ~51.5 days. Between these dips, there 
are periods that appear to be free of extinction lasting 
a few days ea ch, indicative that there may be large gaps 
in the disk. iGalan et al.l ()2010D proposed that similar 
dips in the 2008/09 EE Cep eclipse were due to gaps in 
a multi-ring disk. 

There may also be substructure on timescales shorter 
than a day with variations up to 1 mag on timescales 
shorter than a day. Figures [5] and [5] show detailed struc- 
ture at the beginning and end of the eclipse. If this sub- 
structure is due to the occulting body then it contains 
a remarkable wealth of structure. We have examined 
SuperWASP light curves of neighboring stars of similar 
magnitude, as well as J1407 outside of eclipse, and seen 
no such variations, implying that the hourly variations 
are due to the eclipse. We explore hypotheses to explain 
the eclipses in Sec. 3.6, and develop a model in Sec. 3.7. 

3.4. Rotation Period and X-ray Emission 

Young stars show variability on the timescale of days, 
induced by the presence of starspots on the rotating sur- 
face. Variability at the ^^0.1 mag level can be seen in 
the data, so we carried out a search for periodicity in 
the star's light curve to determine the star's rotational 
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Fig. 8. — SuperWASP {dark filled circles) and ASAS {open tri- 
angles) photometry during the portion of the eclipse immediately 
before the deepest minimum. During some nights the SuperWASP 
photometry shows (unphysical) jumps between two photometric 
levels at the tenths of magnitu de level - a systemat ic effect seen in 
other SuperWASP studies (e.g. [Norton et al.ll20Tll 'l. 
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Fig. 9. — SuperWASP and ASAS photometry for the 2nd part 
of the eclipse; same as for Fig. [S] 

period. The ASAS photometry and SuperWASP pho- 
tometry are measured in slightly different bandpasses, 
and we measure this magnitude difference by taking the 
median magnitude of the data in each data set over the 
2008 season, where there is no sign of long term trends 
in the photometric light curve. We measure a systematic 
offset of 0.143 magnitudes between ASAS and SWASP 
y-band photometry (likely due to the SWASP photom- 
etry being calibrated to the Tycho Vt system), so we 
add an offset to the SWASP photometry to put it on 
the ASAS V system. For the SWASP data, we calculate 
the median magnitude of each night and use this for the 
subsequent analysis. 

We take the photometry from SWASP and ASAS and 
perform a Lomb-Scargle periodogram analysis on both 
data sets. False alarm probabili ties (FAPs) are esti mated 
using the method described in IPress et al.l (|1992[ ) . The 
photometry over the 2008 season is shown in the top 
panel of Figure [TUl along with the Lomb-Scargle peri- 
odograms of the SWASP and ASAS data in the lower 
panels. Both light curves show a highly significant peri- 
odicity of 3.20 days, with FAPs of 10"^ and 10"^ for the 
ASAS and SWASP datasets respectively, with no other 
detectable periods seen over the sampled period ranges. 

The star has an X-ray counterpart in the ROSAT 
All-Sky Survey Faint Source Catalog (RASS-FSC; 
iVoges et al.|[2000l ). with marginally detected flux of fx 
= 0.0359 ±0.0148 ct s'^, and hardness ratioes of HRl 
= -0.04 ±0.42 and HR2 = 0.06 ±0 62. U sing the energy 
conversion factor of lFleming et al.l ()1995D , this translates 
to an X-ray flux in the ROSAT band of fx = 2.9 x 10"" 
erg cm~^, and using our previous distance and bolo- 
metric luminosity estimates, Lx — 10^^'* ergs"^ and 
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\og{Lx/ Lboi) — -3.4 dex. A K5-type star with rotation 
period of 3.20 days would be predicted to have soft X- 
ray emission around the saturatio n level {\og{Lx/ Ltoi) 
~ -3.2 ±0.3: iPizzolato et all 120031 ) . perfectly consistent 
with the observed ROSAT X-ray flux (log(Lx/-^fco/) - 
-3.2), and consistent with other Sco-Cen pre-MS stars 
(|Mamaiek et al.ll2002[ ). 




4S50 4600 
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the results of our analysis. In the upper panel we show 
that we have one or more photometric points for all trial 
periods up to 850 days, and the lower panel shows the 
mean magnitude of the photometric points at those peri- 
ods where we have one or more photometric points. We 
conclude that there is no evidence for any eclipse events 
with photometric coverage up to 2330 days (6.4 years), 
and that there are no eclipse events in any periods up 
to 850 days (2.3 years). Approximately half of periods 
between 850 days and 2330 days are ruled out (Fig. [TT]). 

If the shallow depression seen in early 2001 (see Fig- 
ure [7]) is an eclipse then the period is 6.12 years and just 
about at the 2330 day limit. EE Cep has a similar period 
and duration and a lso exhibits variations in eclipse depth 
(|Galan et al.ll2010D . Hence we should consider the possi- 
bility that J1407 also exhibits large variations in eclipse 
duration and depth. 
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Fig. 10. — Top: Photometry of the star over the 2008 season, 
where the SuperWASP (SWASP) points are daily median values. 
Middle: Lomb Scargle pcriodogram for SuperWASP photometry, 
and Bottom: Lomb-Seargle periodogram for the ASAS photometry. 
False alarm probabilities are indicated with horizontal dashed lines. 
A period of 3.20 days, presumably due to starspots and stellar 
rotation, is detected strongly and independently in both datasets. 
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Fig. 11. — Photometric coverage for periodic events at different 
trial periods. The upper panel shows the periods where we have 
one or more photometric points within an eclipse event, and the 
lower panel shows the mean magnitude in those particular cases. 



3.5. Constraints on the Eclipse Period 

We detect a single deep eclipse in our data set, but 
there is the possibility that another deep eclipse occurred 
during a period when there was no photometric coverage. 
We determine if there could be other eclipses that are 
undetected in the data set by choosing a trial period 
for the eclipses and plotting the light curve modulo this 
period. We count how many folded photometric points 
lie within the deepest part of the known eclipse, and 
determine the mean magnitude and standard deviation 
of these points. 

We define a deep eclipse event as one where the stellar 
magnitude fades by one magnitude or greater. For the 
known eclipse, we estimate that this lasts for 15 days, 
and we set the trial periods P for values starting at 200 
to 2500 days in steps of one day. In Figure [TT] we show 



3.6. Other Explanations 

The case for the primary being a pre-MS K star at 
d ~ 130 pc seems to be secure. The primary exhibits: 
(1) rapid rotation (P = 3.2 days), (2) complimentary 
saturated X-ray emission consistent with the rapid ro- 
tation (typical for young dwarf or pre-MS stars), (3) 
strong Li consistent with other pre-MS Sco-Cen mem- 
bers, (4) proper motion statistically consistent with Sco- 
Cen membership, and predicted kinematic distance har- 
monious with pre-MS status, and lastly (5) spectral ap- 
pearance consistent with being dwarf or pre-MS, but not 
a giant. 

Besides the apparent spectral evidence, there are rea- 
sons to exclude J1407 as a possible giant or supergiant. 
With these observations of the primary in mind, we 
briefly present and pass judgement on several hypothe- 
ses regarding the agent responsible for J1407's unusual 
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eclipses. If J1407 were a K5 giant with absolute mag- 
nitude similar to the K5III standards N Vel and 7 Dra 
(My ^ -l-2£3, then its apparent V magnitude would 
be consistent with a distance of ?a4.3 kpc. J1407's total 
proper motion (32 masyr^^) would then imply a tan- 
gential velocity of ^670 kms^"'^, i.e. faster than the local 
Galactic escape velocity. Photometrically, there is no 
hint of long-term periodicity characteristic of red giants 
(i.e. Mira variability). The situation is worse if the star 
were a K5 supergiant. If it shared the absolute magni- 
tude of the K5 lb standard a CMa (My =i -4.3), then 
J1407 would represent a young, massive star wl8 kpc 
away, ~6 kpc above the Galactic disk midplane, with 
tangential velocity of w2800 kms~^. The presence of 
strong Li absorption. X-ray emission, and a 3-day peri- 
odicity would also seem extraordinarily unusual for a K5 
giant or subgiant. Hence, we rule out J1407 being an 
evolved, giant or supergiant late-K star. 

We discuss some of the possible explanations for the 
observed eclipse, and pass judgement on their plausibil- 
ity. 

• Eclipses by stellar or suhstellar companion alone: 
No plausible stellar or substellar companion can be 
responsible for dimming the K5 pre-MS star J1407 
by more than 3 magnitudes (>95% dimming), and 
eclipses by such an object would not explain the 
irregular shape, depth, and duration of the eclipse. 

• Is the "primary" a red giant that is eclipsing a 
fainter, bluer star? The rare cases of eclipsing bi- 
naries that eclipse by >few magnitudes are usu- 
ally cases of a red giant transiting a smaller, hotter 
dwarf star ( e.g. RV Aps; A2V-fK4III; depth ~ 
1.5 mag; K haliullin et a l. 2006) or symbiotic bina- 
ries (e.g. AR Pav; depth ~ 6 mag, period ~ 604 
day; iQuiroga et al I I2002D . There is no hint from 
the spectrum of J1407 that it could contain a giant 
star, a hot component, or constitute a symbiotic 
binary. As stated before, we think we can safely 
exclude the hypothesis that the J1407 primary is 
an evolved late-K giant or supergiant. Such a sce- 
nario would also not explain the eclipse structure 
in the weeks before and after the main deep eclipse. 

• Could the obscuration be associated with a disk 
orbiting a compact stellar remnant? The system 
is too young to contain a neutron star or white 
dwarf. Given the age of the system (~16 Myr), 
any black hole would have ha d a progenitor mass of 
>14 M© HBertelh et al.ll2009D and would have been 
an extremely large red supergiant and/or Wolf- 
Rayet star before its supernovae. A black hole 
would likely be a m uch stronger s ource of X-rays 
(Lx > 10^2 erg s~^ : IVerbuntl 119931 ) than observed 
(Lx — 10^° ergs~^) if it accreted from a disk. Also, 
the system's proper motion appears to be comov- 
ing with Sco-Cen within ^1 kms~^, so if there were 

Using Hipparcos V magnitudes (Ferryman & ESA 1997) and 
revised Hipparcos parallaxes (van Loeuwcn 2007), we estimate that 
the high qu ality K5 III spectral standards N Vel and 7 Dra 
(IKeenan &: M cNeil 198£) have absolute magnitudes of My — -1.19 
and -1.14, respectively, and the K5 lb stand ard a CMa has Mv — 
-4.3 (assuming extinction ~ 0.12 mag : IBobvlev et al.i r2006). 



a companion that supernovaed and removed a sub- 
stantial amount of mass from the system, then why 
is J1407 not a runaway star? Hence it seems very 
implausible that the obscuration is associated with 
a disk orbiting any type of stellar remnant. 

• Can a circumbinary or circumstellar disk about the 
star explain the obscuration? The pre-MS binary 
system KH 15D has exhibited photometric varia- 
tions and eclipses over the past half century that 
are attrib uted to the effect of a precessing circumbi- 
nary disk (jHerbst et al.ll2010[ ). The eclipses are 3.5 
mag deep in all bands and last for a significant 
fraction (one third or 16 days) of the orbital pe- 
riod (48 days). The single star V718 Persei in the 
young cluster IC 348 also exhibits prolonged (^1 
year) eclipses of about 1 mag with a period of 4.7 
yr. Its eclipses are attributed to an inner edg e of a 
circumstehar disk (lGriiiin_eiaJj |200l |200g|) . For 
both of these systems the eclipses last a signifi- 
cant fraction of the orbital period implying that 
the eclipsing object nearly fills the orbital plane. 
V718 Per has a weak IR excess cor responding to a 
"thin, low-mass disk" (|Grinin et al.l[2 008). KH 15D 
has negligible mid-IR excess (C. Hamilton-Drager, 
priv. comm.). J1407 lacks a near- or mid-infrared 
excess that would indicate a warm dust disk of sub- 
stantial optical depth (and the lack of strong emis- 
sion lines in the spectrum also indicates no evidence 
for an accretion disk). Including the gradual dim- 
ming phase, the eclipse on J1407 lasted about 54 
days. As discussed in the previous section the pe- 
riod analysis suggests that the orbital period P > 
850 days (2.33 yr) so the ratio of eclipse duration to 
orbital period must be less than 0.06. This is signif- 
icantly lower than the ratios for either V718 Persei 
or KH 15D (0.2 & 0.3, respectively) and suggests 
that a circumbinary or circumstellar disk about the 
K star cannot account for the eclipse. 

• Could the eclipse be due to a circumstellar disk that 
occulted the star once due to the relative motions 
of the Sun and J 14-07? We have not yet positively 
identified more than a single eclipse, so at present 
it is possible that the eclipse was a one-time occur- 
rence. As a UCL member, the tangential motion 
of J1407 on the sky is 32 masyr"^ or 20 kms^^ at 
its predicted distance of 128 pc. What if we inter- 
pret the obscuration as being due to a geometrically 
thin circumstellar disk orbiting J1407, with the op- 
tical depth changing due to the motion of the Sun 
relative to J1407? The eclipse depth would then 
suggest an optically thick mid-plane (Amag > 3.5 
mag, or T > 3.2). In a week, the Sun-J1407 line 
only sweeps 0.6 milliarcsecond due to their relative 
motion. Assuming the disk to be of similar size 
to typical planetary orbits ('^10 AU), this would 
translate to sweeping the disk in the z-direction 
approximately ^4(rrfisfe/10 AU) km per week. A 
two week eclipse would correspond to a disk '^10 
km thick, suggestive of a remarkably thin disk with 
aspect ratio (height over radius) of ^10^^, simi- 
lar to the rings of Saturn. Besides the thinness of 
the disk, there are other problems with this see- 
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nario: (1) it does not explain the nearly symmetric 
dimmings at ±12 and ±26 days from the inferred 
eclipse minimum, (2) the similarity with the peri- 
odic eclipsing object EE Cep would have to be co- 
incidental, and (3) the future detection of another 
similar eclipse would obviously negate the idea. 

• Could the eclipses be due to a circumstellar disk 
orbiting a star more massive than the K5 star? 
We can consider the possibility that the system is 
like e Aurigae with the more massive object ob- 
scured by the eclipsing disk. Hiding a dwarf star 
more massive than the K star seen would require 
an e dge-on disk (such as s een in images of HK 
Tau; iStapelfeldt et al.lll998[ ) but again there is no 
evidence for an infrared excess from J1407. The 
12/im WISE flux corresponds to a flux XFx 1.5 x 
IQ-i^W m~^. We have used the 12/xm flux be- 
cause the signal to noise is significantly higher than 
that at 22^m. The error at 12/im is only 3% of the 
flux so at best a disk could be emitting with a 12/xm 
infrared flux of XFx 5 x 10"^'' W m"^. If only 
1/10*'* of the disk luminosity is emitted at 12/im 
then the total infrared flux of our source is at most 
Fir < 5 X 10"i^Wm"^. The solar luminosity 
at a distance of 128 pc (that estimated for J1407) 
corresponds to 2 x 10"-^^ W m~^. Thus the total 
infrared luminosity is at most Ljr < 2.5 x 
It would be difficult to hide a main sequence star 
and remain below this luminosity, typical of debris 
disks, even if the disk were an edge-on transition 
disk with an inner hole. It is more likely that the 
object hosting the occulting disk has a lower mass 
than the K star. In this case the disk infrared lu- 
minosity could be consistent with this upper limit. 

3.7. A Preliminary Model 

We now consider the possibility that the eclipse could 
be due to occultation by a circumsecondary or circum- 
planetary dust disk with the secondary object and its 
disk in orbit about the K5 primary star (analogous to 
what has been proposed for EE Cep). The period anal- 
ysis above suggests that the orbital period P > 850 day 
(2.33 yr), which for nii = 0.9 M© suggests an orbital ra- 
dius of >1.7 AU and an orbital velocity of <21.7 kms~^. 
The hypothetical disk appears to produce some obscura- 
tion during a minimum of At « 54 days, and hence the 
obscuration occurs during a faction fecUpse <6.3% of the 
period. Using equation [S] for the fraction of time spent 
in eclipse we find that this limit implies that 



using an eclipse time of 54 days and period of 12 years) 



TO2 < 21MjC 



(8) 



where we have used mi = 0.9Mq estimated for J1407. 
If the disk fills the Hill or tidal radius (C ~ 1) then the 
secondary is likely to be a brown dwarf. The secondary 
could have a higher mass if the disk radius only partly 
fills the tidal radius estimated from its semi-major axis. 
If the secondary is in an eccentric orbit and its disk is 
truncated tidally at pericenter then ^ estimated from the 
semi-major axis would be lower than 1. 

If the dimming seen in 2001 corresponds to a secondary 
eclipse then the fraction of the period spent in eclipse is 
f eclipse ~ 0.024 Corresponding to (using equation [5] and 



TO2 - l-2Mj^- 



(9) 



For ^ ^ 0.2 expected for circumplanetary disks this gives 
m2 ~ O.IM0 and an M dwarf. Both this mass estimate 
and the previous one estimated from the 850 day period 
limit suggest that the mass of the companion must be low 
and in the brown dwarf or low mass M star regime. The 
longer the period, the lower fecUpse, and the lower the 
estimated mass of the secondary. However if the period 
is longer then the separation between primary and sec- 
ondary would be larger making it easier to resolve using 
a high angular resolution imaging system. For example, 
a separation of 10 AU at a distance of 130 pc corresponds 
to 77 mas, resolvable perhaps with ape rture-masking in- 
terfer ometry on large telescopes (e.g. [Ireland &i KrausI 
[2OOI) . 

The period and length of the eclipse can be used to 
estimate the disk radius. Inverting equation [S] 
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As expected, only large objects are capable of causing 
such a long eclipse. The disk radius is only weakly de- 
pendent on the period (proportional to P~^/'^) and would 
be smaller for a more distant companion. 

Regions during the eclipse with little dimming could 
be int erpreted in terms of gaps in the disk (as by 
iGalan et al. 2010 for EE Cep). The gaps (regions la- 
belled Zi, Z2, Yi and Y2 on Figure [5]) each last a few 
days. The ratio of the width of these gaps to one half of 
the eclipse time is approximately 3/26 or 0.11 suggesting 
that the disk must have an aspect ratio h/r smaller than 
this ratio. If we assume that the gaps are twice the Hill 
radius of an objected embedded in the disk then the ratio 
of the gap to total eclipse time gives an upper limit on 
the ratio of the third power of the ratio of the satellites 
to secondary mass. 



< 



eclipse 
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where uis is the mass of the gap opening satellite. For 
tgap/teciipse ~ 0.06 this corresponds to ms/'rn2 ^ 10"'^. 
If the secondary has a mass of IMj then the satellites in 
the disk could have mass lower than that of Earth, and if 
the secondary is an M star of mass O.IMq then the gap 
opening satellites could have mass lower than Saturn. 

Photometric variations on daily timescales (e.g., see 
Figures [3 [9]) suggest that there are hourly variations in 
the eclipse depth. The ratio of a few hours (or a quarter 
day) to the half echpse length is about 0.01 suggesting 
that the disk has an extremely low aspect ratio of h/r < 
0.01. If so then the disk could not be a gaseous disk but 
must be a planetesimal disk or a ring system. 

We are currently attempting to model the eclipse light 
curve in terms of an optically thick inner disk which 
caused the primary deep eclipse and a system of rings 
of lower optical depth which caused the smaller dips in 
the weeks before and after the primary eclipse (Scott, 
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Moolekamp, & Mamajek, in prep.). We have had suc- 
cess modeling the primary echpse, however the periph- 
eral dips have been more difRcult to model. Light from 
the primary star is modeled as a collection of spherically 
symmetric points (iWilson fc DevinnevI 119711 ) and a sec- 
ondary star or planet modeled as a sphere (we assume 
the tidal deformations of both objects to be negligible 
because the length of the eclipse and minimum period 
implies a large disk and distance to the primary ). Lim b 
darkening was calculated following ' van Hammd (pM ). 
The disk and rings are assumed to be thin debris disks 
of dust with uniform density and opacity. The input pa- 
rameters were the mass of the primary, the orientation 
of the debris disk, the orientation, period and eccentric- 
ity of the orbit, the inner and outer radius of each ring, 
and the opacity of each ring. Physical parameters as 
a fu nction of age and m ass for the primary were taken 
from lBaraffe et all (Il998f) an d for the low-mass compan- 
ion from lBaraffe et all (|200"^ . 

A preliminary good fit (but by no means unique) to 
the daily-averaged light curve for J1407 is shown in Fig. 
1121 s-nd the model parameters are listed in the caption 
and summarized in Table [H A diagram showing the ge- 
ometry of this preliminary good fit is shown in Fig. 1131 
Besides a thick inner disk needed to explain the deep 
eclipse labelled C in Fig. 5, our toy model includes two 
"rings" of different optical depth for explaining features 
Bi and B2 in the light curve. An additional optically 
thin outer ring is needed to explain features Ai and A2 
in the light curve. This outer ring is not included in this 
model, but has the following approximate parameters: 
Tin ~ 200 Rc (where Rc is the companion radius, for this 
model assumed to be 1.46 Rj„p ~ 104,000 km). Tout — 
250 Rc, r_L ~ 0.09. Another gap is needed to explain 
maxima Zi and Z2 in Fig. 5, with inner and outer radii 
of approximately ~163 Rc and ~200 Rc. 

Approximately how much dust could be responsible for 
the eclipses that we are seeing? The densest of Saturn's 
named rings - the B ring - has a radially-averaged optical 
depth of orde r unit y and surface density of ~ 50 gcm~^ 
(iZebker et al.lll985t ' lde"Pater fc Lissaueill2001l : rTiscarend 
|2012| ). implying an approximate opacity of k ~ 0.02 
cm^g~^. Adopting this mass opacity for J1407's dust 
annuli, the thick inner disk would have a mass of ~0.8 
Mj\/oon, the two rings shown would have masses of ~0.2 
and '-~^0.1 Mj\,foon, and the outermost ring would contain 
~0.5 Mmoou of dust mass. Our toy model assumes that 
the companion and its disk system is situated 8.7 AU 
from J1407, however this is only definitively constrained 
to be >1.7 AU given the time-series photometry avail- 
able. So any estimates of the physical scale and mass of 
the disk system will obviously scale with the companion's 
orbital separation from J1407. For k — 0.02 cm^ g~^, the 
total ring mass ranges from ^3.6 Mjv/oon for P = 2.33 yr, 
to '--^O.S Mmoou for P — 200 yr. Over this same range, 
the outer edge of the outermost ring scales from 60 mil- 
lion km (0.4 AU) for P = 2.33 yr, down to 14 million 
km (0.09 AU) for P = 200 yr. We are in the process 
of improving the code, attempting to fit the intra-night 
light curves, and making predictions of brightness of the 
companion's disk in the infrared, in order to give further 
constraints on the size and orientation of the debris sys- 
tem, and the physical parameters of the gaps and dust 
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Fig. 12. — A simple, non-unique model attempting to fit the gross 
features of the nightly mean SuperWASP photometry for J1407. 
This model contains around an object orbiting J1407 girded by a 
thick inner disk, a gap, and two rings with a small gap between 
them. The thick inner disk is used to model the deep eclipse feature 
"C" in Fig. 5, and the two rings are modeled to fit the features 
"Bi" and "B2" in the same figure. The companion object has a 
test orbital inclination of 89°. 955, axial tilt with respect to the 
orbital plane of 13°, orbital period 9862 days (a = 8.7 AU), and 
radius Rc = 1.46 Rj^p. The model contains a thick inner disk with 
Tj_ = 0.5 and outer radius 76 Rc, a first "ring" with optical depth 
Tj_ = 0.2 between 106 and 127 Rc, and a second "ring" with optical 
depth Tj_ = 0.05 between 128 and 163 Rc. Yet another outer ring 
is needed to fit dips before and after the time range plotted. 




Fig. 13. — Diagram of J1407's dust disk model whose correspond- 
ing light curve is plotted in Fig. 1121 The K5 star J1407 is plotted 
to scale on the right (radius = 0.96 R0). The thick (r = 3) inner 
disk is needed to produce feature "C" in Fig. 5, and the tilt is 
largely responsible for the asymmetric eclipse. The two annuli of 
lower optical depth (rperp = 0.2 and 0.05) are used to fit features 
Bi and B2 in the light curve. The gap between the thick inner 
disk and inner ring is needed to model features Yi and Y2 in Fig. 
5. The outer ring hypothesized to explain features Ai and A2 in 
the light curve in Fig. 5 is not included in this diagram. The 
intra-night SuperWASP light curves are suggestive of much more 
substructure within the "rings" than represented here. 



4. THE PROBABILITY OF SEEING ECLIPSES BY 
CIRCUMSECONDARY OR CIRCUMPLANETARY DISKS 
IN A SAMPLE OF YOUNG STARS 

We first estimate the probability of a circumplanetary 
disk eclipse using probability distributions for giant plan- 
ets estimated from radial velocity surveys. We then es- 
timate the probability of a circumsecondary disk eclipse 
based on surveys of young binary stars. 

4.1. Probability of seeing eclipses by circumplanetary 

disks 

To estimate the probability of detecting a circumplan- 
etary disk eclipse we must consider the number of stars 
that host gas giant planets. The period and mass distri- 
bution of gas giants estimated from Doppler radial ve- 
locity surveys is 



dN 



CM-3.1±0.2p0.26±0.1^lQg^^^lQgp 



(12) 
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TABLE 5 

Ring Model Parameters 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Ring 




Rout 


T 


T± 


Mass 




(Rc) 


(Rc) 






Mjv/oon 


"Rochester" 


1? 


81 


3.0 


0.5 


0.8 


"Sutherland" 


106 


127 


1.0 


0.2 


0.2 


"Campanas" 


128 


163 


0.3 


0.05 


0.1 


"Tololo" 


200 


255 


0.5 


0.09 


0.5 



Note. — Ri„ is inner radius, Rout is outer radius, 
T is optical depth through ring along line of sight 
between observer and the primary star, is the 
optical depth through the ring perpendicular to the 
ring place. Ring sizes are parameterized by compan- 
ion radius, assumed to be 1.46 Rjup (104,000 km). 
This model assumes a companion orbital period of 
27 years (a = 9.7 AU) and ring opacity of k = 0.02 
cm^ g~ ^ . To scale the radii for different assumed or- 
bital periods, multiply them by factor (P/27 yr)~^/'^ 
(where P > 2.33 yr). To scale the ring masses, mul- 
tiply by (P/27 yr)-2/3(K/0.02 cm^g-l)-!. Ring 
nicknames come from locations where observations 
were taken or analysis carried out for this study 
(The SuperWASP observations were taken using the 
SuperWASP-South observatory located at Suther- 
land, South Africa. The ASAS survey was carried 
out at Las Campanas Observatory, Chile. Spectra 
of the host star were taken with the SMARTS 1.5- 
m on Cerro Tololo. Discovery and analysis of the 
system took place at University of Rochester) . 

where the normahzation constant C is such that the frac- 
tion of FGK stars with a planet in the mass range 0.3 - 
10 Mj (where Mj i s a Jupiter mass) and p eriod range 2- 
2000 days is 10.5% (|Cumming et al.ll2008r ). In units cor- 
responding to measuring planet masses in Jupiter masses 
and orbital periods in days, the value of the normaliza- 
tion is C = 1.4 X 10"'^. Integrating over the masses we 
find dN = C"P°-26d log P with C = 4.5 x IQ-^. The dis- 
tribution gives a probability of a FGK star hosting a gas 
giant between 1-5 AU (365-1825 days), interior to this (2 
to 365 days), and from 5-20 AU, in each case of about 
fg ~ 0.05 or 5%. 

We consider a sample of stars restricted so that they 
have already depleted circumstellar disks (i.e. are post- 
accretion pre-MS stars) but are young enough that 
they could host circumplanetary disks of sufficient op- 
tical depth to produce detectable eclipses. One could 
choose a sample based on eliminating stars with ev- 
idence of accretion and selecting for age based on 
chromosp heric activity, cluster or association member- 
ship (e.g. IMamaiek et al]|2002l : IMamaiek fc Hillenbrandl 
l2008f ). Typical subgroups of OB associations have '-^lO^ 
stars, and ages of ~ 3-20 Myr (e.g. dc Zccuw et al. 1999; 
IBriceno et al.l 120071. during which the majority of stars 
have just rec ently ceased ac creting from circumstel- 
lar disks (e.g. iMamaiefl |2009() . Light curves of post- 
accretion ~1 Mq pre-MS stars in the nearest OB as- 
sociations (e.g. Sco-Cen, Ori OBI, etc.) can be searched 
in existing SuperWASP and ASAS datasets, and indeed 
such an effort is currently underway by our group. 

If each star in the sample is observed a single time the 
fraction that would be observed in eclipse we estimate 
as /i f orient f eclipse fg whcre wc multiply the number 
of systems with gas giants, fg, by the fraction of orbit 
spent in eclipse, fecUpse (equation [5]), and fraction of ori- 



entations capable of giving eclipse f orient (equation [3]). 
For our three ranges of semi-major axis radii 

f,^e^^'"i-'"^''yfs (13) 

^lo-mv^yf— Vf^) f M 

V0.2y \mj) \Mq) Vo.5/vo-05y 

This probability is very low and implies that high- 
cadence monitoring of many stars is required to detect 
circumplanetary disk eclipses. 

We now consider the same sample of stars but con- 
tinuously monitor them throughout the planet's orbital 
period. The fraction (for each range of semi-major axis) 
that would exhibit an eclipse of a circumplanetary disk 

would be /c ~ f orient fg, 

Restricting our study to planets between 1-5 AU (with 
fg ^ 0.05) and monitoring them for 10 years (approxi- 
mately the orbital time at 5 AU for ~1 Mq) the above 
fraction fc suggests that if we monitor 10"* stars 10^ years 
old then ^2 of them should exhibit circumplanetary disk 
eclipses. 

A system with a recurring circumplanetary eclipse 
would make it possible to study eclipses in depth, so 
discovery of systems with short orbital periods are im- 
portant. For example, circumplanetary disk substructure 
could be inverse modeled by observed high cadence light 
curves during eclipse. Unfortunately due to the smaller 
Hill radius size closer to the star, we expect that the life- 
time would be short for circumplanetary disk in a smaller 
orbit about the star. One could search for circumplane- 
tary disk eclipses in systems that have not yet lost their 
outer circumstellar disks (for example systems such as 
/3 Pic or DM Tau). However a circumstellar disk is a 
large object and the probability that a circumplanetary 
disk occults the star but the circumstellar disk does not 
occult the star is probably even more negligible. 

4.2. Probability of Seeing Eclipses by Circumsecondary 

Disks 

We consider a 10 year photometric survey of a sam- 
ple of weak-lined T Tauri stars. The fraction that would 
exhibit an eclipse by a disk would be the fraction that 
are binaries (~ 0.5; iDuchene et"all 120071: IKraus et all 
1201 If ) times the fraction that have binary periods less 
than 10 years, times the fraction that have a secondary 
with an optically thick disk which could produce a sig- 
nificant eclipse, times the probability that such systems 
are oriented in a way giving eclipses (given by equation 

ED. 

For young binary systems in star-forming regions the 
estimated fraction of mixed systems with primary a 
weak lined T Tauri star and the secondary a classical 
T Tauri star is not low, and could be as large as ^ 1/3 
onin et afl 12007^ . The fraction of young binary sys- 
tems that have a weak lined T Tauri primary and a sec- 
ondary with a pass ive non-accre t ing d isk is lower; the 
survey described by iMonin et al.l ()2007D contains only 1 
out of about 80 binaries. A tidally truncated disk around 
a low mass secondary is expected to ha ve a shorter ac- 
cretion lifetime than the primary's disk (|Armitage et al.l 
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I1999D . though its planet formation timescale could be 
longer. Thus a high mass ratio binary (with mass ratio 
of order g ~ 0.1) with a classical T Tauri phase primary 
and a passive disk about the secondary should be rela- 
tively rare. A mid-IR survey of about 65 binary young 
stars finds that abou t 10% contain passive dust disks 
(|McCabe et all [2006h and for one of these the disk is 
ho sted by the second ary. Thus of the binaries studied 
bv lMonin et"all pOOTf ) and lMcCabe et all (|2006[) we can 
crudely estimate that 1/100 could be like J1407 with a 
weak-line T Tauri primary and a low mass secondary 
with a passive disk (however this fraction should dimin- 
ish as one gets to older pre-MS stars with a g es of >10'' 
yr) . As described by iPrato fc Weinbergeil (|2010( ) , bi- 
nary systems with primaries "...classified as weak-lined 
T Tauris, unresolved, might also harbor truncated disks 
around the secondary stars. Such small structures could 
go undetected as the result of dilution from a relatively 
bright primary. Circumstellar disks with central holes 
that show excesses in the mid-infrared but not in the 
near-infrared, and which do not show signatures of ac- 
cretion, may be present but are effectively undetectable." 
J 1407 may be in this class and an example of a relatively 
rare young binary with weak-lined T Tauri primary and 
low mass secondary hosting a passive disk. 

The number of b inaries is flat in l og pe r iod or semi- 
major axis space (|Halbwachs et al.l 120031 : iKraus et al.) 
[20T1 and young bin aries are similar to field popula- 
tion in this respect (jPuchene et al.l [20071) . Based on 
this distribution about 1/3 of all binaries have periods 
less than 10 years (dividing the semi-major axis ranges 
into three ranges 1-10 yr, 10-100 yr, and 100-1000 yr). 
Thus the number of weak-lined T Tauri stars that are 
binaries with periods less than 10 years and contain low 
mass secondaries with passive disks would be of order 
1/600. Using a mass ratio of 0.1, y — 0.3 and ^ = 1 
we estimate the fraction oriented such that they can 
give eclipses forient ^ 0.1. Altogether a 10 year sur- 
vey of weak-lined T Tauri stars might have a probability 
of detecting an eclipse of order 1/6000 giving a similar 
probability to that estimated above for circumplanetary 
disk eclipses. An estimate folding the mass distribution 
and lifetime distributions would improve this estimate. 
However this is difficult to formulate as binary identifi- 
cations are not com plete at mass ratios less than 0.1 (e.g., 
IKraus et al.ll201l[l and the number of secondaries host- 
ing passive disks in the 10^ year old age range is not well 
characterized. There are hints that the protoplanetary 
disk fraction decay timescale is systematically longer for 
lower-mass stars and br own dwarfs com pared to Sun-likc 
stars and massive stars ()Mamaiek " 2009) , however further 
observations will be useful in constraining these findings 
for components of binary systems. 

5. CONCLUSION 

In this study we have estimated the probability that 
a system hosting a gaseous circumplanetary or circum- 
secondary disk about a planet could occult a star. The 
existence of circumplanetary disks after the dissipation 
of the protosolar (circumstellar) nebular disk has been 
postulate d from formation scenarios for the Galilean 
satelUtes (iCanup fc \yardll2002l: iMagni fc Coradinill2004l: 
iWard fc Canud |2010[ ). Because such a disk would be 
large the probability that a system hosting one is ori- 



ented in such a way that it can occult the star is tiny, 
but not zero. Because the lifetime of a circumplane- 
tary disk could be longer for planets at large semi-major 
axis, light detected from outer exoplanets such Fomal- 
haut b may arise from such a disk. We estimate that 
eclipses from the thick inner circumplanetary disks that 
spawn regular satellite systems around gas giants may 
last for days, however tenuous outer disks of lower op- 
tical depths to larger fractions of the Hill radius could 
persist for weeks, depending upon the planet's mass and 
semi-major axis. We estimate that a survey monitor- 
ing 10'' stars that are approximately 10^ years old for 10 
years would likely yield at least a few circumplanetary 
and circumsecondary disk eclipse candidat es. The 8.4- 
m La rge Synoptic Survey Telescope (LSST; llvezic et"al] 
l2008f l will photometrically monitor many thousands of 
pre-MS stars in catalogued OB associations, young star 
clusters, and star-forming regions during its scans of the 
Galactic plane region. The LSST survey design should 
enable estimation of sub-mas/yr proper motions which 
will allow kinematic membership assignments of newly 
discovered pre-MS stars to young clusters and associa- 
tions of determinable distance and age. Optical spec- 
troscopic follow-up will be necessary to confirmation of 
youth (via Li absorption, c.f. our study of J1407) and 
estimation of reddening for confirming that the star has 
luminosity and isochronal age consistent with the other 
cluster/association members. Even given LSST's pro- 
posed field revisit time (^3 days, but twice per night), 
a J1407-like eclipse would have been easily detected. 
Higher cadence follow-up imaging from smaller dedicated 
telescopes will be necessary for detailed characterization 
of the eclipse light curves, but LSST monitoring of young 
stellar groups should yield some candidate disk eclipse 
objects. 

In a survey of a few hundred 10^ year old stars we have 
discovered a deep long eclipse in 2007 on the pre-main 
sequence K dwarf star J1407. This star was selected to 
be in the appropriate age range and could host either 
a circumplanetary disk or a lower mass secondary star 
with a disk. The lack of infrared emission suggests that 
the mass of the unseen object is much lower than that of 
the solar-mass K star. Limits for the period and the ob- 
served eclipse time suggest that the unseen object host- 
ing the disk is low mass, perhaps in the substellar regime. 
Substructure in the eclipse suggests that the disk is thin 
h/r < 0.01 and has gaps that may contain satellites with 
mass ~10~'^ times that of the secondary. Follow up with 
a radial velocity study is important as radial velocity 
measurements could put limits both on the period and 
mass of the secondary. 

The complex eclipse of J1407 that took place in 2007 
is slightly asymmetric and contains significant substruc- 
ture, similar to the eclipses of the Be star EE Cep that 
have be en interpreted in te rms of an occulting planetary 
system (jGalan et al.ll20ldl ). As has been proposed for 
this system, the asymmetry of the eclipse could be due 
to the impact paramet er of the d isk with respect to the 
line of sight (.Mikolaiewski fc Graczvk.l999i) . Variations 
in eclipse depth in this system are attributed to a pos- 
sible third comp anion that tilts t he orbital plane of the 
eclipsing system (jTorres fc Stefan ik 2000). Identifying a 
second eclipse for J 1407 will allow measurement of the 
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eclipse period, and planning of observing campaigns sim- 
ilar to those launched for EE Cep and e Aurigae. 

Constraints from the gas giant satellite systems in our 
own solar system suggest that their circumplanetary disk 
structures could have produced quite complex eclipses 
if seen in transit, with dense inner regions, gaps where 
satellite formation is taking place, and low density disks 
possibly extending to large fractions of the Hill radius. 
Such eclipses seen among young stars may provide re- 
markable laboratories for testing satellite and planet for- 
mation scenarios. Regardless of the nature of the disked 
companion of J1407 (low-mass star, brown dwarf, or 
gas giant planet), detailed observations of future eclipses 
should provide useful constraints on either circumsec- 
ondary or circumplanetary disk structure, and the early 
evolution of planets and/or satellites. 

We have used data from the WASP public archive in 
this research. The WASP consortium comprises of the 
University of Cambridge, Keele University, University 
of Leicester, The Open University, The Queen's Uni- 



versity Belfast, St. Andrews University and the Isaac 
Newton Group. Funding for WASP comes from the 
consortium universities and from the UK's Science and 
Technology Facilities Council. The star exhibiting the 
unusual eclipse was discovered in a spectroscopic sur- 
vey using the SMARTS 1.5-m telescope, and the sur- 
vey and support for EM and MP were funded by NSF 
award AST-1008908. EM, MP, FM, and ES acknowl- 
edge support from the University of Rochester College 
of Arts and Sciences. AQ acknowledges support through 
NSF award AST-0907841. This research has made use of 
the NASA/IPAC Infrared Science Archive, which is op- 
erated by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the National 
Aeronautics and Space Administration. EM also thanks 
David James, John Subasavage, Andrei Tokovinin, War- 
ren Brown, and Catrina Hamilton-Drager for discussions, 
and Fred Walter for scheduling queue observations of the 
star and standards on the SMARTS 1.5-m. 

Facilities: SuperWASP, ASAS, CTI0:1.5m, 
CTI0:1.3m 
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